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THREE-DIMENSIONAL PERSPECTIVE SOFTWARE 
FOR REPRESENTATION OF DIGITAL IMAGERY DATA 

By Bobby G. Junkin'* 

SUMMARY 

A generalised three-dimensional perspective software 
capability has been developed at the NASA/NSTL Earth Resources 
Laboratory within the framework of a low-cost computer-oriented 
geographically based information system using the Earth Resources 
Laboratory Applications Software (ELAS) operating subsystem. 

This perspective software capability was developed primarily 
to support data display requirements at the NASA/NSTL Earth 
Resources Laboratory, It provides a means of displaying three- 
dimensional feature space object data in two-dimensional picture 
plane coordinates and makes it possible to overlay different 
types of information on perspective drawings to better understand 
the relationship of physical features. An example topographic 
data base is constructed and is used as the basic input to the 
plotting module, Examples are shown which illustrate oblique 
viewing angles that convey spatial concepts and relationships 
represented by the topographic data planes. 


*Earth Resources Laboratory, National Space Technology 
Laboratories, National Aeronautics and Space Administration, 
NSTL Station, Mississippi 39529 
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THREE-DIMENSIONAL PERSPECTIVE SOFTWARE 
FOR REPRESENTATION OF DIGITAL IMAGERY DATA 


INTRODUCTION 

The NASA/NSTL Earth Resources Laboratory (ERL) has 
developed several computer software programs and procedures 
(references 1-4) within the framework of a computer-oriented 
information system (references 5-6) for the processing and 
analysis of data from disparate, geographically oriented 
base maps and from remote sensor aircraft and satellite 
systems. These different types of information are compiled 
into data bases which contain information on land use, 
elevation, slope, soil series, rainfall, population density, 
etc. The capability to manipulate, store, analyze, display, 
and disseminate the large volumes of data in these data 
bases has evolved through research efforts at the ERL 
(References 5 and 6) , 

Examples of several practical modular systems, with 
emphasis on low cost, are given in reference 5. These 
systems, an example of which is shown in figure 1, consist 
of an image display system, a graphic digitizer, a small 
digital computer, and an output recording device. All 
hardware components used in these low-cost data processing 
systems are off-the-shelf. The software consists of a 
Landsat multispectral scanner data reformatting program, 
a series of supervised and unsupervised spectral-pattern- 
recognition programs, a program to reference the image 








data to a map base, a data storage and retrieval program, 
and various applications programs. 

The software support concept has evolved into an operating 
subsystem referred to as the Earth Resources Laboratory Appli- 
cations Software (ELAS) (reference 7) , This software system 
accepts as input a variety of data types including topographic 
data tapes from the National Cartographic Information Center 
(NCIC) (reference 8) . Users of these data require oblique 
viewing angles to convey spatial concepts and relationships 
represented by the topographic data planes. For example, the 
relationship of Landsat classification components to elevation 
may be visualized by perspective color mapping of the classifi- 
cation data plane on the topographic elevation data plane. 

Other graphic information such as aspect or slope can be 
placed on the perspective plot to relate the data elements 
to the terrain profile. This report addresses the procedures 
involved in the development of a generalized three-dimensional 
perspective software capability to support the above display 
requirements and presents typical results using a topographic 
data base . 
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THREE-DIMENSIONAL GRAPHICS WITH PERSPECTIVE 


General Overview 


There are several steps required to determine the two- 
dimensional perspective picture plane coordinates for the 
display of three-dimensional feature space object data 
(Reference 9). First, the coordinates of the feature space 
object are translated to the vantage point coordinate system. 
Then the vantage point coordinates are rotated through two 
angles to align with the vantage point line of sight passing 
through the origin of the reference coordinate system. Finally, 
the perspective transformation is applied to obtain the two- 
dimensional picture plane coordinates for the actual display 
on two-dimensional display devices (i.e., X-Y plotters on 

CRT screens) . An integral part of this display is the removal 
of hidden lines on the other side of the object from the 
vantage point. 


Coordinate Systems 

There are three rectangular coordinate systems involved 
in the derivation of two-dimensional picture plane coordinates. 
These are the reference coordinate system (X,Y,Z) , the vantage 
point coordinate system (X‘ ,Y',Z’), and the picture plane 

coordinate system (X M ,Y" ,Z M ) . A brief description of each 
follows . 
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1. (X,Y,Z) Reference System . This is a right-handed fixed 
system of axes. The coordinates of the vertices of all 
objects as well as the vantage point (X 0 ,Y 0 ,Z 0 ) are given 
in this system. 

2. (X'tY'.Z*) Vantage Point System . The vantage point 
(X 0 ,Y o ,Z o ) specified in the reference system is taken as 
the origin of the vantage point coordinate system, in 
which coordinates are denoted as (X'jY’.Z 1 ). This is a 
left-handed coordinate system wherein the Z’ axis goes 
from (X 0 ,Y 0 ,Z 0 ) through the origin of the (X,Y,Z) reference 
system. The X* axis is parallel to the XY plane of the 
(X,Y,Z) reference system. 

3. (X^Y 11 ) Picture Plane System . This is a two-dimensional 
coordinate system which represents the plane of the actual 
perspective drawing. This plane goes through the origin 
(X,Y,Z) of the reference system and is chosen perpendicular 
to the Z* axis. The X" axis is chosen to lie in the XY 
plane. These coordinate systems are depicted in Figure 2. 

Coordinate Transformations 

We first consider the translation from the reference 
coordinates to the vantage point coordinates . From Figure 3 
we see that: 
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Figure 2. REFERENCE, VAN * AGE POINT AND PICTURE PLANE COORDINATE SYSTEMS 
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Figure 3. INITIAL POSITION OF REFERENCE AND VANTAGE POINT 
COORDINATE SYSTEMS 
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These equations can be expressed more conveniently in matrix 


form as : 


m »m 

X' 


-1 

0 

0 

Y' 


0 

0 

1 

Z' 


0 

-1 

0 


or: 

X = % X 



(3) 


Two rotations are now required so that the vantage point 
coordinate system is oriented whereby the Z r axis goes through 
the origin of the reference coordinate system. This is 
accomplished with a rotation about the Y' axis by the angle a, 
followed by a rotation about the X* axis by the angle 3. The 
rotations about the X 1 and Y' axes are shown in Figures 4 and 5 , 
respectively. The final rotation in Figure 5 results in the Z' 
axis passing through the origin of the reference coordinate 
system. A rotation about the Y* axis through the angle a results 
in the following transformation matrix: 


cos G 

0 

m 

-sina 

0 

1 

0 

sinu 

0 

cosa 


(4) 


After this rotation the Z' axis points at (0,0,Z). 

The second rotation operates on the result of the first 
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Figure 4. ROTATION THROUGH ANGLE a ABOUT Y' AXIS 
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rotation and results in the Z' axis going through (0,0,0), 
This transformation matrix is*. 



10 0 
0 cos# sin3 

0 -sin3 cos3 

tm 


Thus, the total transformation can be written as: 


X = R-£ Ry R-p X 


( 5 ) 


( 6 ) 


Expanding the above equation yields : 


■» 

X' 


-(X-X 0 ) cosot + (Y-Y 0 )sina 

Y' 

cr 

- (X-X 0 )sin3sinot - (Y-Y 0 ) sin3cosa 4* (Z-Z Q )cos3 

Z' 


- (X-X 0 ) cosPsinot - (Y-Y 0 ) cosPcosa - (Z-Z 0 )sin3 


From basic trigonometry, the following relations are obtained 
from figure 4 : 


sina » X 0 /d 
cosa = Y 0 /d 
and from figure 5 : 

sin8 = Z q /D 
cos3 ~ d/D 


( 8 ) 


(9) 
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where : 
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( 10 ) 


Projective Transformation 


Once the Z’ axis is oriented to pass through the origin 
of the reference system, the 3-D coordinates (X'jY’.Z') must 
be projected onto the X"Y"-picture plane to generate the 
actual 2-D coordinates. This projection is shown in figure 6. 
Noting similar triangles in this figure, the point V n (X l * l Y M ) 
in the picture plane is related to the point V f (X 1 , Y 1 ,Z ' ) in 
the vantage point system by the following equations: 


X M = DX’/Z’ 
Y M » DY’/Z’ 


(ID 


Scaling 

The final step before displaying the (X",Y") perspective 
data on a CRT screen or other display device is to scale the 
data to fit within display limits. Basically, this involves 
scanning the (X",Y") data to find maximum and minimum X" and Y H 
values for use in computing linear transformation of the form 


12 



i. 

OF POOR QUALITY 



/ 


Figure 6. 


GEOMETRY FOR PERSPECTIVE PROJECTION 




t *. V 

OF POOR 


n 



$ “ A o v " + h- For elevation data obtained from contour 
profiles we want A 0 and to be such that Y^ n is 0 and Y” 
is some value Z max , Thus, we have: 


A o Y min + A X « 0 
A o Y max * A 1 “ Z max 


( 12 ) 


Solving for A 0 and A^: 

A o " ^max/ ( Y max “ Y min^ 
A 1 " A o Y min 


(13) 


The elevation data values are then computed from: 

[WC^ax - ^min>] ( Y " - Y min) < 14) 


Hidden Line Procedure 

There are several algorithms in the literature for 
eliminating hidden lines from the perspective representation 
of data in three-space (references 10-14) , The two-dimensional 
representation of a three-dimensional surface consists of 
line segments of a succession of curves, The basic approach 
used herein is to eliminate all of those line segments which 
are behind other surfaces and which would not be visible in 
a plot of the data in two dimensions. A lower boundary (or 
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horizon array) for all lines is set up below which no line 
is drawn. Visibility arrays are set up and then the line 
segment points of the curve closest to the observer are 
tested. If the point is visible then it is inserted in the 
visibility array and becomes a point in the horizon array. 

If it is not visible then the next point of the line segment 
is tested and the process repeated until all points along 
the curve have been tested. The visibility array now contains 
the points for plotting the first curve in the foreground. 

This logic is then applied to the line segment points of the 
next curve furthest from the observer. At each step the 
horizon array is also tested and updated. 
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INITIAL APPROACH 


Background 

Most large-scale computer facilities have some form of 
software capability for the graphic display of three-dimensional 
feature space data in two-dimensional form. Initial investi- 
gations at the NASA/ERL revealed that there was no 3-D software 
capability available from local vendors for implementation on 
the NASA/ERL computer system. Subsequent efforts led to the 
development and implementation of a basic 3-D program for a 
specific function. 

Program Implementation 

The complete program as developed on the NASA/ERL com- 
puter system for a specific function is shown in figure 7. 

Buffer arrays are defined as shown in this figure. The 
DATAPLOT software routines PLOTS, PLOT and DPINIT are used to 
create the plot file for plotting on the STATOS plotting system. 

An offset origin on the STATOS is defined and boundary vectors 

and loop counters are initialized. The logical parameters A, B, 
and S(I,J) are used in the "hidden-line" logic whereby a visi- 
bility check is performed on each point along a curve to determine 
if it is hidden by other parts of the surface. The domain of the 
function Z - F(X,Y) is also defined as shown in this figure. 

Results 

One example of a surface that can be represented in a 
perspective display plot is defined by the following 
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Figure 7. FLOW CHART FOR THREE-DIMENSIONAL PERSPECTIVE SOFTWARE USING A 
SPECIFIC FUNCTION 










































single-valued continuous exponential function: 


Z 


EXP — 225 (2X—Y) 


3.239 C'Y- 


•. 5X) 


( 15 ) 


An example of a perspective function that can be used to 
generate (X',Y') data for 2.-D display plotting is given by: 


X” = Z' + X’ cosa 
Y" = Y' + X' sina 


(16) 


An example plot using these functions is shown in figure 8. 

This plot was created using a function domain of -6 £ (X* ,Y’)£ 6, 
This example serves to illustrate one approach for the display 
of three-dimensional data in two-dimensional perspective. An 
obvious shortcoming to this approach is the function definition 
requirement. A much more flexible and efficient approach is 
described in the next section. 
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THREE-DIMENSIONAL SOFTWARE DEVELOPMENT 


Disc ussion 

The development of three-dimensional perspective software 
capability to support the display requirements addressed 
earlier in this report has been accomplished under the ELAS 
operating subsystem. A functional diagram of this geographically 
based information system is shown in figure 9. The generalized 
three-dimensional program as developed for running on the NASA/ERL 
computer system is summaiized in the flow chart shown in figure 
10. Various functional programs such as CFSUB, RDWR, RIO, 

ILBYTE, and ISBYTE were used in this development. 

Description of Program 

The application module name for this program is referred 
to as PDDD. This application module determines the two- 
dimensional perspective picture coordinates for the display 
of three-dimensional data, such as the NCIC elevation data 
that has been processed through the topographic application 
modules (reference 15) . This program provides for the overlay 
of other types of data onto this perspective. The program 
can be executed either in the demand (interactive) or batch 
mode of operation. The data are processed beginning with the 
surface profile closest to the observer. Each profile is 
read and scaled according to a previously determined scale 
factor. Overlay data are then extracted from the appropriate 
input. The (X",Y") picture plane data and the overlay data 
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Figure S. GEORAPHICALLY BASED INFORMATION SYSTEM 
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Figure 10. FLOW CHART FOR GENERALIZED THREE-DIMENSIONAL PERSPECTIVE SOFTWARE 
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are then passed to the subroutine HIDLIN (figure 1L) to 
determine if any portion of the line is visible, Visible 
segments of the line are then passed to the subroutine INTERP 
(figure 12) which writes the data to a disk file for subsequent 
plotting on an X-Y type plotter or viewing on an image display 
device. 


Directives and Operating Instructions 

The input for PDDD is a standard ELAS data file referred 
to as ID1, The output is another ELAS data file known as ODF. 
Two channels of data from ID1 are specified as input channels; 
these two channels may or may not be distinct depending on 
what type of output the user wants. For example, suppose 
the user has a two- channel input file where channel 1 contains 
NCIC elevation data as output from T0P6 . Then each value in 
channel 1 represents a range of elevations. Suppose channel 
2 contains the Landsat classification values that correspond 
to the elevation values for each pixel in channel 1. If the 
user specifies the input channels as 1 and 2, the program would 
output the classification value contained in channel 2 at the 
3-D perspective coordinate (line, element and elevation value) 
determined by the corresponding data of channel 1. If the 
user specifies the input channels as 1 and 1, the program 
would output the elevation value in channel 1 at the 3-D 
perspective coordinate determined by the same data (line, 
element, and elevation value) in channel 1. 
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The following directives control the operation of the 
program: 

LD - List directives. 

SP - Set parameters, 

PH - Print header information (input and output) . 

PF - Prepare file. Write zeroes in output channel. 

PD - Plot data. (PD, R, T where R * rotation angle 
in degrees and T = tilt angle in degrees.) 

The following parameters, with typical values in parentheses, 
are required input: 


IL 

- 

Initial line (1) . 



LL 

- 

Last line (<IL -f 1024), 



IE 

- 

Initial element (1) , 



LE 

- 

Last element (<IE + 1024), 



NCH 

- 

Number of input channels (2) , 



CH 

- 

Channel numbers for input (1,2). 



XNCL 

- 

Input line increment (1) , 



INCE 

- 

Input element inci'ement (1) , 



XOF 

- 

X offset to plot origin (10) . 



YOF 

- 

Y offset to plot origin (10) . 



NGO 

- 

Number of output channels (1) . 



CHO 

- 

Channel number for output (1) . 



SD1S 

- 

Maximum scaled data value (50) . 



X 

- 

Number of rasters per element on 

X-axis 

(2). 

Y 

- 

Number of rasters per element on 

Y-axis 

(4). 
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TOPOGRAPHIC DATA BASE CONSTRUCTION FOR 3-D INPUT 


Digitized topographic elevation data are obtained from 
the National Cartographic Information Center (NCIC) in the 
planar format. These data are digitized from 1 : 250 ,000-scale 
maps and are available from NCIC on 800- or 1600-BPI tapes. 

These data serve as input to a series of three topographic 
application modules for the construction of a topographic data 
base and subsequent input to the three-dimensional plotting 
module PDDD. The first module consists of six overlays --T0P0, 
TOPI, T0P2, TOP 3 , TOP 4 , and T0P5; a second module (T0P6) 
computes slope, aspect, and slope length on the output of 
T0P5; a third module (T6CH) computes the average north-south 
slope and the average east-west slope in addition to the four 
variables computed in T0P6 . Basically, these modules reformat 
the NCIC data from tape into the ELAS data file format, rotate 
the data 90 degrees to a north-south orientation, compute 
mapping coefficients that relate transposed plate coordinates 
to UTM coordinates (eastings and northings) , resample the 
16-bit data to a UTM grid file, and use this file to compute 
slope, aspect, slope length, and average slopes. These 
modules are summarized in figure 13. A technical writeup 
of these procedures with a topographic data base example is 
contained in reference 15 . 
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Figure 13. NASA/ERL TOPOGRAPHIC DATA PROCESSING AND ANALYSIS SOFTWARE OVERVIEW 
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RESUL TS USING OLYMPIC NATI ONAL PARR TO POGRAPHIC DATA BASE 

The TOPQ application modules were, used to construct a 
multi -channel topographic data base Tile containing elevation, 
slope, aspect and slope length parameters for the Olympic 
National Park in the State of Washington at a cell si&e of 
50 meters * 50 meters* Another channel in this data base 
contains the results of classifying a Lands at MSS data set 
for this area and subsequently improved by using elevation 
data. The park encompasses some one million acres and includes 
the Olympic Mountains, Elevations in the Olympic Park range 
from vsea level to approximately 8,000 feet with variations 
in vegetation responding to severe changes in terrain and 
rainfall patterns within the park, 

A vegetation- type land cover classification derived from 
the Lands at’ MSS data and improved using NCIG elevation data 
is shown in figure 14, Further details on the derivation of 
these results are given by Oihula in reference 16, The color 
codes for the 21 classes are shown in figure 15. The relation- 
ship of the classification components to the elevation com- 
ponents may he visualised from figure 16 which contains the 
perspective color mapping of the classification data plane 
on the topographic elevation data plane. This perspective 
plot was obtained using a viewing rotation angle of -10 degrees 
and a tilt angle away from the viewer of 45 degrees. A 
backside view of a section of the elevation and classification 
data pianos for this same area is shown in figure 17, 




ORIGINAL PAS^ 

color photograph 


F ig jre 15. COLOR CODES FOR OLYMPIC NATIONAL PARK VEGETATION CLASSIFICATION DERIVED 

FROM LANDSAT MSS DATA 


CLASSIFICATION COLOR 
COOI 


DESCRIPTION 


1 

2 
3 


5 

6 

7 

8 

9 

10 

II 

17 

13 

14 

15 

16 

17 

18 

19 

20 
21 







WATER 

HARDWOODS OR HARDWOOD/CONIFEROUS MIX 

LOWLAND HARDWOOD OR HARDWOOD/CONIFEROUS MIX; MAT INCLUDE 
BRUSH OR SLASH 

AGRICULTURAL FIELDS WITHOUT VEGETATION OR W'TH 

VERY LITTLE (24); CLEAR CUT AREAS; SLASH; CLEAR CUTS 
WITH EARLY REGENERATION (< a 20 30% VEGETATION COVER 
OR LESS). CONCRETE. OTHER INERTS (INCLUDES SOME URBAN) 

PASTURE A OTHER GRASSLANDS 

BEACH. RIVER GRAVELS. WET SAND DRIFTWOOD. MODERATELY 
DARK COLORED INERT MATERIALS 
SHADOWS 

SNOW ICE OR GLACIERS 

SCRUB HIGH ELEVATION HARDWOOD; VINE MAPLE A SITKA ALDER 
ALPINE MEADOWS (CLASS 44 WITH SPARSE VEGETATION) 
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Figure 17.- A backside perspective view. 



Figure 18 shows a section of this same area with no rotation 
and 60 degrees tilt. The displays in figures 16, 17, and 18 
can be located in figure 14 by the indicated boundary markings. 
These results are but a few examples of how three-dimensional 
perspective graphics can be used to convey spatial concepts 
and relationships represented by topographic data planes. 


34 



35 


Figure 18.- Perspective relationship of classification to elevation for no rotation and 
a tilt of 60 degrees. 



CONCLUDING REMARKS 


A generalized three- dimensional perspective software 
capability has been developed at the NASA/ERL within the 
framework of the Earth Resources Laboratory Applications 
Software (ELAS) operating subsystem, This perspective 
software capability provides a means of displaying three- 
dimensional feature space object data in two-dimensional 
perspective picture plane coordinates and makes it possible 
to superimpose disparate or spatial information on perspective 
drawings to better understand relationships of physical features. 

Digitized topographic elevation data obtained from the 
National Cartographic Information Center are used as input 
to a series of topographic application modules for the 
construction of a topographic data base. These data base 
parameters are subsequently input to the three-dimensional 
plotting module and serve to illustrate one type of data 
base input. The relationship of a Landsat classification to 
elevation is illustrated by the perspective color mapping of 
the classification data plane on the topographic elevation 
data plane. The example displays illustrate oblique viewing 
angles that convey spatial concepts and relationships repre- 
sented by the topographic data planes. 
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